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Abstract 

This note is devoted to the explicit construction of a functional defined 
on all pairs of L 1 functions with small total variation, which is equiv- 
alent to the L 1 distance and non increasing along the trajectories of 
a given system of conservation laws. Two different constructions are 
provided, yielding an extension of the original stability functional by 
Bressan, Liu and Yang. 
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1 Introduction 

Let the smooth map h W 1 define the strictly hyperbolic system of 
conservation laws 

d tU + d x f(u) = o (i.i) 

where t > 0, x € K and u € f2, with f2 C W 1 being an open set. 

Most functional theoretic methods fail to tackle these equations, es- 
sentially due to the appearance of shock waves. Since 1965, the Glimm 
functional [13] has been a major tool in any existence proof for 111. IB and 
related equations. More recently, an analogous role in the proofs of contin- 
uous dependence has been played by the stability functional $ introduced 
in HOI EI] , see also [I] . The functional $ has been widely used to prove the 
L 1 -Lipschitz dependence of solutions to (and related problems) from 
initial data having small total variation, see for example [TJ [2j [TUl [HI 13 
116] . Special cases comprising data with large total variation are considered 
in [HI [HI [T71 [181 CEH]- Nevertheless, the use of <3? is hindered by the neces- 
sity of introducing specific approximate solutions, namely the ones based 
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either on Glimm scheme |13j or on the wave front tracking algorithm [U [12] . 
The present paper makes the use of the stability functional $ independent 
from any kind of approximate solutions. The present construction allows 
to simplify several parts of the cited papers, where the presentation of the 
stability functional needs to be preceded by the introduction of all the ma- 
chinery related to Glimm's scheme or wave front tracking approximations, 
see for instance [9]. 

We extend the stability functional to all L 1 functions with sufficiently 
small total variation. This construction is achieved in two different ways. 
First, we use general piecewise constant functions and a limiting procedure, 
without resorting to any sort of approximate solutions. Secondly, we exploit 
the wave measures, see [U Chapter 10] and give an equivalent definition that 
does not require any limiting procedure. Furthermore, we prove its lower 
semicontinuity. 

With reference to [I] for the basic definitions related to (jl.ip . we state 
the main result of this paper. 

Theorem 1.1 Let f generate a Standard Riemann Semigroup S on the 
domain T>g defined at (|2.6p . Then, the functional 3 defined at (|3.4p by 
means of piecewise constant functions coincides with the one defined at (I4.5P 
by means of wave measures. Moreover, it enjoys the following properties: 

(i) 3 is equivalent to the L 1 distance, i.e. there exists a C > such that 
for all u,u £ T>s 



(ii) 3 is non increasing along the semigroup trajectories, i.e. for all u,u G 
T>s and for all t > 



(Hi) 3 is lower semicontinuous with respect to the L 1 norm. 

Taking advantage of the machinery presented below, we also extend the 
classical Glimm functionals [H [13] to general L 1 functions with small total 
variation and prove their lower semicontinuity, recovering some of the results 
in [3], but with a shorter proof. 

2 Notation and Preliminary Results 

Our reference for the basic definitions related to systems of conservation 
laws is [3]. We assume throughout that G £1, with Q open, and that 




u — it|| L i < 3(u, u) < C ■ \\u — u 



L 1 • 



S(S t u,S t u) < S(u,u) . 
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(F) / G C 4 (f2;M n ), the system (II. ip is strictly hyperbolic with each char- 
acteristic field either genuinely nonlinear or linearly degenerate. 

Let Xi(u), . . . , X n (u) be the n real distinct eigenvalues of Df(u), indexed so 
that Xj(u) < Xj + i(u) for all j and u. The j-th right eigenvector, normalized 
as in (|2.ip - (|2.2p . is rj(u). Let a i— > Rj(a)(u), respectively cr i— > Sj(a)(u), be 
the rarefaction curve, respectively the shock curve, exiting u, so that 

dRj(a)(u) , , w .> , 
V =rj(A»H) • (2.1) 

If the j-th field is linearly degenerate, then the parameter a above is the 
arc-length. In the genuinely nonlinear case, see [U Definition 5.2], we choose 
a so that 

r)\ ■ r)\ ■ 

^ ( Rj (a)(u)) = k, , -^L(S 3 (o){u))=k 3 (2.2) 

where k±, . . . ,k n are arbitrary positive fixed numbers. In [1] the choice kj = 1 
for all j = 1, . . . , n was used, while in (2' another choice was made to cope 
with diagonal dominant sources. The choice (12,2ft preserves the properties 
underlined in [4J Remark 5.4] so that the estimates in [H Chapter 8] still 
hold. Introduce the j-Lax curve 



a i ^ ipj (a) (u) 



Rj(cr)(u) if cr > 
Sj{a)(u) if a < 



and for cr = (ci, . . . , er n ), define the map 

*(<r)(w _ ) = Vn(o"n) o . . . o ^(ctx)^ - ) . 

By [H § 5.3], given any two states u~,u + G 17 sufficiently close to 0, there 
exists a map i£ such that 

(<7i, . . . ,cr„) = E(u~,u + ) if and only if u + = ^{a)(u~) . (2.3) 

Similarly, let q = (qi, . . . , q n ) and define the map S by 

S(q)(«") = 5 n ( 9n ) o . . . o Sxiq^iu') (2.4) 

as the gluing of the Rankine-Hugoniot curves. For any two states u~ ,u + as 
above, there exists a unique q such that u + = S(q)(u~). 

Let u be piecewise constant with finitely many jumps and assume that 
TV(u) is sufficiently small. Call T(u) the finite set of points where u has 
a jump. Let a x ^ be the strength of the i-th. wave in the solution of the 
Riemann problem for (|l.ip with data u(x-) and u(x+), i.e. (<r x ,i> • • • > &x,n) = 
E (u(x— ), u(x+)) . Obviously if x ^ 1{u) then a Xj i = 0, for all i = 1, . . . , n. 
As usual, ^4. (it) denotes the set of approaching waves in it: 



((x,i),(y,j)) G (Z(n)x{l,...,n}) 2 : 
x < y and either i > j or i = j, the i-th field 
is genuinely non linear, min {cr^j, a y j } < 
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As in [13J or [4, formula (7.99)], the linear and the interaction potential are 



V(u) = XT X^l "^! and Q( n ) = XT \ <T x,i (T V,3\ 

x£l(u) i=l ((x,i),(y,j))eA(u) 



Moreover, let 

T(n) = V(u) + C • Q(u) (2.5) 

where Co > is the constant appearing in the functional of the wave-front 
tracking algorithm, see [U Proposition 7.1]. Recall that Co depends only on 
the flow / and on the upper bound of the total variation of initial data. 

Remark 2.1 For fixed x\ < . . . < xn+i, the maps 

(«!,..., U N ) I ^ Qp = i!i«X[x [tA+1 [ / 

are Lipschitz continuous. Moreover, the Lipschitz constant of the maps 

is bounded uniformly in N , a and u a for a ^ a. 
Finally, for 5 > sufficiently small, we define 

T>1 = |w € L 1 (R, U) : v is piecewise constant and Y(v) < <5 j 

V s = dV*s 

where the closure is in the strong L 1 -topology. Unless otherwise stated, we 
always consider the right continuous representatives of maps in T>s and T>%. 
For later use, for u G T>s and r\ > 0, introduce the set 

B v (u) = |« € L 1 (R; Q): u G Dj and ||v - u|| Ll < 77} . (2.7) 

Note that, by the definition of T>$, B v (u) is not empty and if r/i < 772, then 
Bm(u) Q B m (u). Recall the following fundamental result, proved in [BJ: 

Theorem 2.2 Let f satisfy (F). Then, there exists a positive 5 Q such that 
the equation (jl.ip generates for all 5 G ]0, S [ a Standard Riemann Semigroup 
(SRS) S: [0, +oo[ x T>$ 1— > which is Lipschitz in L 1 . 

We refer to [U Chapters 7 and 8] for the proof of the above result as well 
as for the definition and further properties of the SRS. 
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3 The Piecewise Constant Functions Approach 

Extend the Glimm functionals to all iiGPjas follows: 

Q(u) = lim inf Q(v) and T (u) = lim inf Y(v) . (3.1) 

r)-+0+ v eB v (u) V-*0+veB v (u) 

The maps r\ — > inf^ gB r v \ Q(v) and r\ — > inf^g^ ^ T(v) are non increasing. 
Thus the limits above exist and 

Q(n) = sup inf Q(v) and Y(u) = sup inf T(t>) . 

ri>0v£B v (u) n>0 v£B v (u) 

We prove in Proposition 13.41 below that Q, respectively Y, coincides with Q, 
respectively T, when evaluated on piecewise constant functions. Moreover, 
we prove in Corollary 14.41 that Q also coincides with the functional intended 
in (5j formula (1.15)] or [U formula (10.10)]. Preliminarily, we exploit the 
formulation (|3.ip to prove directly the lower semicontinuity of Q and T. 

Proposition 3.1 The functionals Q and Y are lower semicontinuous with 
respect to the L 1 norm. 

Proof. We prove the lower semicontinuity of Y, the case of Q is analogous. 

Fix u in D$. Let u v be a sequence in T>s converging to u in L 1 . Define 
£v = \\u u — "u|| L i + l/v. Fix v v £ B ei/ (u u ) so that 

T(w„)< inf T(v) +e u < T(u u ) + e u . 

v€B £v (u v ) 

From \\v v — u|| L i < \\v v — u u \\ L i + \\u u — u\\ L i < 2e u we deduce v u £ B2 £u {u) 
and the proof is completed with the following estimates: 

inf T(v) < T(v v ) < T(u v ) + e v ; 

T(u) = lim inf Y(v) < lim inf T(it I/ ) . 

The next proposition contains in essence the reason why the Glimm 
functionals Q and T decrease. Compute them on a piecewise constant 
function u and "remove" one (or more) of the values attained by u, then 
the values of both Q and T decrease. 

Let u = J2aei Ua X[x a ,x a+1 [ be a piecewise constant function, with u a € 
fi, x\ < X2 < ■ ■ ■ < xn + i and / be a finite set of integers. Then, we say 
that Ui,«2) • • • ) un is the ordered sequence of the values attained by u and 
we denote it by (u a : a G J) . 

Proposition 3.2 Let u,u G T>%. If the ordered sequence of the values at- 
tained by u is (u a : a E I) and the ordered sequence of the values attained by 
u is (u a :a G J), with J C J, then Q(n) < Q(u) and T(u) < T(n). 
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Ua-1 Ua Ua+l 



Figure 1: Proof of Proposition 13.21 u a is attained by u and not by u. 

Proof. Consider the case jjl = ft J + 1> see also [H Lemma 10.2, Step 1]. 
Then, the above inequalities follow from the usual Glimm interaction esti- 
mates p3] , see Figure CD 

The general case follows recursively. □ 

The next lemma is a particular case of [H Theorem 10.1]. However, the 
present construction allows to consider only the case of piecewise constant 
functions, allowing a much simpler proof. 

Lemma 3.3 The junctionals Q and Y ; defined on T>%, are lower semicon- 
tinuous with respect to the L 1 norm. 

Proof. We consider only T, the case of Q being similar. 

Let u v be a sequence in V* s converging in L 1 to u = Y^a u aX[x a ,x a+1 [ G 
as v — > +oo. By possibly passing to a subsequence, we may assume that 
Y(u u ) converges to liminf^ +00 Y(it„) and that u u converges a.e. to u. 
Therefore, for all a = 1, . . . , N, we can select points y a 6 ]x a , x a +i [ so that 
lim u ^ +00 u u (y a ) = u(y a ) = u a - Define 

u u = ^2u u (y a )X[x a ,x a+1 [- 

a 

By Proposition 13.2} Y(u„) < Y(u u ). The convergence u u (y a ) — > u a for all 
a and Remark 12.11 allow to complete the proof. □ 

Proposition 3.4 Let u£T>* s . Then Q(u) = Q(u) and T(u) = Y(u). 

Proof. We consider only Y, the case of Q being similar. 

Since u € T>t, we have that u € B„(u) for all rj > and Y(-u) < 
Y(u). To prove the other inequality, recall that by the definition (|3.ip of 
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Y, there exists a sequence v u of piecewise constant functions in T>% such 
that v v — ► u in L 1 and T(i?^) — > T(u) as v — > +00. By Lemma 13. 3| 
T(u) < liminf^^+oo T(v u ) < Y(-u), completing the proof. □ 

Therefore, in the sequel we write Q for Q and T for Y. 

Since we will need the explicit dependence of the Stability Functional on 
the various quantity it is made of, we introduce the following notations. If 
5 € ]0,<5 o [, for any v £ denote by a x ^ the size of the i-wave in the 
solution of the Riemann Problem with data v(x— ) and v(x+). Then define 

Aj[v](x) = ^2\a yJ \, Af[v](x) = ^2\d- y>j \, for j = l,...,n. 

y<x y>% 

If the i-th characteristic field is linearly degenerate, then define Aj as 

l<j<i i<j<n 

While if the i-th characteristic field is genuinely nonlinear 

l<j<i i<j<n 

+ A tM(x) ■ X P)+00[ (?) + A^M(x) ■ X yoofl[ (q) 
Now choose v, v piecewise constant in T>'* s and define the weights 

Wi[v,v](q,x) = 1 + KiAi[v){q,x) + KiAi[v]{-q,x) 

+/si«a (Q(«) + Q(«)) • 1 ' ' 

the constants k\ and K2 being defined in [U Chapter 8]. We may now define a 
slightly modified version of the stability functional, see [3 [201 EI] and also [U 
Section 8.1]. Namely, we give a similar functional defined on all piecewise 
constant functions and without any reference to both e-approximate front 
tracking solutions and non physical waves. 

Define implicitly the function q(x) = (qi(x), . . . ,q n (x)) by 

v(x) = S (q(x)) (v(x)) 

with S as in (|2.4p . The stability functional $ is 

n f+OO 

&(v,v) = Y \qi(x)\ -Wi[v,v] (qi(x),x) dx. (3.3) 

i=l J -°° 

We stress that <& is slightly different from the functional <!> defined in [H 
formula (8.6)]. Indeed, here all jumps in v or in v are considered. There, 
on the contrary, exploiting the structure of e-approximate front tracking 
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solutions, see [U Definition 7.1], in the definition of <J> the jumps due to 
non physical waves are neglected when defining the weights Ai and are con- 
sidered as belonging to a fictitious (n + l)-th family in the definition [U 
formula (7.54)] of Q. To stress this difference, in the sequel we denote by 
$ £ the stability functional as presented in [U Chapter 8]. 

Remark 3.5 For fixed x\ < . . . < x/v+i 3 x\ < ... < x N , 1; the map 



Ui,..., Up} 
Ui,...,U N 



\ / TV N ^ 

^ * '^2 u aX[x a ,x a+1 [:'^2 u aX[x a ,x c , +1 [ 
/ \a=l a=l i 



is continuous. Indeed, since both maps q i— > (/Xj + ^(q) an d Q *—* QX^ qXq) 



are Lipschitz, for any fixed i£l the integrand in (|3.3p depends continuously 
on {u a }^ = i, {u a }a=i an d the Dominated Convergence Theorem applies. 

We now move towards the extension of to T>$. Define 

& v (u, u) = inf {&(v, v):v G B v (u) and v G ^(u)} 

The map i] — > H r? (u, m) is non increasing. Thus, we may finally define 

H(u, u)= lim H^(n, u) = sup H r? (u, u) . (3-4) 

We are now ready to state the main result of this section. 

Theorem 3.6 Let f satisfy (F). The functional S:T>s x T>s *—> [0,+oo[ 
defined in (|3.4p enjoys the properties (i), (ii) and (Hi) in Theorem \l.l[ 



Here and in what follows, we denote by C positive constants dependent only 
on / and Sq. We split the proof of the above theorem in several steps. 

Lemma 3.7 For all u, u G T>%, one has H(u, u) < *&(u, u). 

Proof. By the definition (|2.7p we have u G B v (u) and u G B v (u) for all 
n > 0, hence S^(u, it) < u) for all positive r?. The lemma is proved 
passing to the limit rj —* 0+. □ 

Lemma 3.8 Let u, u G T>$ and q G M. Assume that u is given by 

N 



U = 

a=l 



where x\ < . . . < xjy + i and y a G [x a ,x a+ i[ are given real points. Then, 
Ai [u] (q,x) + K 2 Q{ii) < Ai [u] (q,x) + K2Q(u) + C ■ \\u(x) - u(x)\\ . 
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Figure 2: Exemplification of point 2 in the proof of Lemma 



Proof. Fix x G R and prove the above inequality passing from u to u 
recursively applying three elementary operations: 

1. it?' is obtained from w only shifting the position of the points of jump 
but without letting any point of jump cross x. More formally, if w = 

EaWaXKa.k+iI Wittl ^ a < W ' = Ea™aX[^ +1 [ With < £a+l 

and moreover x € [£a,£a + i[n [£ a ,£a+i[ f° r a suitable a, then 



Aj [u/] (g, x) + K 2 Q(V) = Aj [w] (q, x) + k 2 Q( 



w 



Indeed, if all the jumps stay unchanged and no shock crosses x, then nothing 
changes in the definition of Aj and Q. 

2. w' is obtained from w removing a value attained by w on an interval 
not containing x, see Figure More formally, if w = w aX[£ a ,£ a +i[ w ith 
£a < (,a+i and x L then 



IV 



w. 



or 



In both cases, 



w 



Aj [to'] (g, x) + re 2 Q(u/) < Aj[u;](g, x) + fc 2 Q(io) . 

Indeed, consider for example the situation in Figure [2 The two jumps at 
the points and in w are substituted by a single jump in w' at the 
point £a+i- The points and are both to the right of x, therefore 
the waves in w' at the point £a+i which appear in Ai[w'](q, x) are of the 
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same families of the waves in w at the points £q and £a+i which appear in 
A.i[w](q,x). Since all the other waves in Aj are left unchanged we have 

n 

Ai[w'](q,x) - Ai[w](q,x) < ^ |^ a+1J - - o^ yj - a ?a+1J 

3=1 

where 

a L+l,j = E j ,W' (Ca+1+)) = Ej (Wa-l,W a+1 ) 

°W = Ej (w(^-) ,w(Ca+)) = Ej (Wa-l,Wa) ■ 

Therefore, the increase in Aj evaluated at x is bounded by the interaction 
potential between the waves at and those at (,a+i and is compensated by 
the decrease in k 2 Q, as in the standard Glimm interaction estimates. 

3. w' is obtained from w changing the value assumed by w in the interval 
containing x. More formally, if w = ^2 a w a X[^ a: ^ a+1 [ with £ a < £, a +i and 
x G [£a,£a+l[> then 

w' = WaX^a+xl + U 4X[f a ,&+i[ • 

In this case 

Ai [w 1 ] (q, x) + k 2 Q{w') < Ai[w](q,x) + k 2 Q{w) + C ■ \\wa - w' a \\ 

< Ai[w](q, x) + k 2 Q(w) + C ■ \\w{x) — . 

This inequality follows from the Lipschitz dependence of Ai[w](q, x)(x) and 
Q(w) on Wa with a Lipschitz constant independent from the number of 
jumps, see Remark 12. II 

Now for x G [xa,Xa+i [ we can pass from u to the function w defined by 

w = Y u(y a )x[ Xa , Xa+1 [ + u(x)x[ x& , Xik+1 [ 

applying the first two steps a certain number of times. We obtain 

Ai[w](q,x) + k 2 Q(w) < Ai[u](q,x) + k 2 Q(u) . 
Finally with the third step we go from w to u obtaining the estimate: 

Ai[u](q, x) + K 2 Q(n) < Ai[w](q,x) + k 2 Q(w) + C ■ \\u(x) — w(x)\\ 

< Ai[u](q,x) + k 2 Q(u) + C ■ \\u{x) - u(x)\\ 

which proves the lemma. □ 
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Lemma 3.9 Let u, u, u, u €T>t. Assume that u and u are given by 



N 



N 



x a+1 [: 



U 



a=l 



a=l 



where xi < . . . < x N+1 , y a € [x a , x a+1 [, xi < . . . < Xfi +1 , y a G [x a ,x a+1 [ 
are fixed real points. Then, 



& (u,u) < $ (u, u) + C I \\u - n|| L i + 



u — u 



L 1 



Proof. Introduce q(x) = (qi(x), . . . , q n (x)) and q(x) = (qi(x),..., q n (x)) by 
u(x) = S (q(x)) {u(x)) and u(x) = S (q(x)) (u(x)) with S defined in (12, 4|) , 



* (^-u, uj — 3? (u, u) 

n r+oo . - 

/ 1 ^(^I^K"] (Qi( x )> x ) - \qi(x)\Wi[u,u] (qi(x),x) \ 
^ / | ft 0*01 {Wj[w,w] (%(a;),x) - Wi[u,u] (qi(x),x) } dx 

._1 J — oo 



da; 



8=1 

+ 



i=l 



Since the map g i— » q- Wj[u, x) is uniformly Lipschitz, the first integral 
is bounded by 



+oo 



\qi(x) - qi(x)\ dx < C I \\u - u\\ L i + 



u — u 



i=l 

Concerning the second integral, observe that by Lemma [37 
Wi[u,u] (qi(x),x) -Wi[u,u] (qi(x),x) 

< k\C < \\u(x) — u(x)\\ + u(x) — u(x) 



L 1 



and, since u?j(x) is uniformly bounded, the Lemma is proved. 



□ 



Proposition 3.10 The functional <1? ; defined on T>% x T>%, is lower semi- 
continuous with respect to the L 1 norm. 

Proof. Fix u, u in V* & . Choose two sequences of piecewise constant maps 
u u ,u u in D| converging to u, u in L 1 . We want to show that $>(u,u) < 
liminf^^+oo <&(u u , u u ). Call I = liminf v _» +00 <&(«„, u v ) and note that, up 
to subsequences, we may assume that lim^+oo <&(u u ,u v ) = I. By possibly 
selecting a further subsequence, we may also assume that both u u and u u 
converge a.e. to u and u. 
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Let {xi, . . . , a^Ar+i} be the set of the jump points in u and u and write 

N N 



U 



Xa -)£a-\-l [ ' 



a=l 



For all a, select y a G ]x a ,2 a +i[ so that as v — > +00, the sequence u u (y a ) 
converges to = u a and u u (y a ) converges to u(y a ) = u a . Introduce 

the piecewise constant function u v = Yl a u v (y a ) X[x a ,x a+ i[- Let u u be de- 
fined analogously. Observe that u v and u v converge pointwise and in L 1 to 
respectively u and u. Remark 13.51 implies lim v ^ +00 $ (u v ,ujj = <&(u,u), 
while Lemma 13.91 implies 



&(u v , u„) < * (u v , u v ) + C - u u \\ Ll + 



L 1 



Therefore, passing to the limit v — > +00, the proof is completed: 

*(u,-u) = «&(«,«) - ^(itj,,^) + ^(itj,,^) - ^(itj,,^) + <&(u u ,u u 
< &(u,u) - &(u u ,u u ) + C ^\\ii u - u,y\\ L1 + 
+&(u v ,u v ) . 

Lemma 3.11 For all u,u € one has H(u, u) > <fr(it, u). 



□ 



Proof. By the definition (|3.4p of S, for all u,u € X>|, there exist sequences 
tV, of piecewise constant functions such that for v — > +00 we have tv — > u, 
v u — ► u in L 1 and <&(uj,, {;„) — > H(it, u). Hence, by Proposition 13.101 



<&(u, u) < liminf Uj/) = H(u, u) , 

V— >+oo 

completing the proof. □ 
Lemma 13.71 and Lemma [3.111 together yield the following proposition. 

Proposition 3.12 For all u, u in T>$, one has S(u,u) = <&(u,u). 

Thanks to the definition (13. 4p of H, we obtain 

Proposition 3.13 The functional S: T>$ x T>s 1— > M is lower semicontinuous 
with respect to the L 1 norm. 
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Proof. Fix u and u in T>s- Let u u , respectively u v , be a sequence in V$ 



converging to u, respectively u. Define e v 



u\ 



+ \\u u 



u\\ 



+ l/v. 



Then, for each v, there exist piecewise constant v v G B £v {u L/ ), respectively 
v v G B £i/ (u u ), such that 



®(v v ,v u ) < S £u (u v ,u v ) + e u < S(u u ,u u ) + e v 



(3.5) 



Moreover 



IL 1 

II 1 



< 
< 



U u\\j_,l 4" 1 1 lit 

^^IIl 1 4" 11^ 



u|| L i < 2e u 
tt|| L i < 2e v 



so that G B 2£v (u) and Uj, G B 2£v (u). Hence, "3 2Ev {u,u) < &(v u ,v u ). 
Using (j3.5fl . we obtain H2 £l/ (u, u) < S(u u ,u u ) + e^. Finally, passing to the 
lower limit for v —* +oo, we have &(u,u) < liminf„^ +00 S(u u ,u u ). □ 



In the next proposition, we compare the functional $ defined in (|3.3[) 
with the stability functional <I> £ as defined in [H formula (8.6)] 

Proposition 3.14 Let 5 > 0. Then, there exists a positive C such that for 
all e > sufficiently small and for all e- approximate front tracking solutions 
w(t,x),w(t,x) o/jnj 



* (w(t, -),w(t, •)) - $ £ (w,w)(t) <Cs- \\w(t, •) - w(t, -)|| L1 • 



Proof. Setting w(t,x) = S (q(t, a;)) (w(t, x)) and omitting the explicit time 
dependence in the integrand, we have: 



*(w(t, •),*(*, 0) -$ e K*)(i) 



< 



i=l 



Wi[ro,w] (%(x),x) -Wi(x] 



dx . 



We are thus lead to estimate 

Wi[w,w] (qi(x),x) - Wi(x 

< ki Ai[w] (qi(x),x) + Ai[w] (-qi(x),x) - Ai(x) 
+kik 2 \Q(w) - Q{w)\ + kik 2 \Q{w) - Q(w)\ . 

The second and third terms in the right hand side are each bounded as in [U 
formula (7.100)] by C e. Concerning the former one, recall that, except when 
qi(x) = or on a finite number of points where w or w have jumps, A^ and 
Aj differ only in the absence of non physical waves in Ai. In other words, 
physical jumps are counted in the same way in both Ai and Aj while non 
physical waves appear in Aj but not in A{. Therefore, the former term is 
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almost everywhere bounded, when qi{x) 7^ 0, by the sum of the strengths 
of all non physical waves, i.e. C e by [4, formula (7.11)]. Finally, using [H 
formula (8.5)]: 



C 



\w(t, x) — w(t, x) II < 



1=1 



i(t, x) I < C ■ \\w(t, X) — w(t, x) I 



(3.6) 



we complete the proof with the following estimate: 

„ n 

$>(w(t,-),w(t,-)) -$ £ (w,w)(t) <Ce V \qi(x)\ dx < C e\\w 

Jm i=1 



w 



L 1 ' 



□ 



Proof of Theorem ] 3. (k The estimates [H formula (8.5)] show that $ is equiv- 
alent to the L 1 distance between functions in Z?|. Indeed, if 5 is sufficiently 
small, then Wj G [1,2] for all i = 1, . . . , n and all x G R, so that 



1 

c ' l|v 



u|| L i < *(u,u) < 2C 



m — v\ 



(3.7) 



C 
1 



To prove (i), fix u, fx G and choose t> £ B v (u), v G B v (u). By (13.71) . 

kx — fx || L i-2?7) < *(«,«) < 2C- (||«-fi|| L i +277) 
( , J«-«|| L i -2r?) < 3 v (u,u) < 2C ■ (\\u - u\\ Ll + 2 V ) . 

The proof of (i) is completed passing to the limit r/ — > 0+. 

To prove (ii), fix u, u € £>,5 and rj > 0. Correspondingly, choose G 
B n (u) and f)^ G B^u) satisfying 

H(u, fx) > S^(u, fx) > $(^, v v ) - rj . (3.8) 

Let now e > and introduce the e-approximate solutions v £ and u£ with 
initial data v £ (0,-) = v ri and uf (0, •) = v v . Note that for e sufficiently small 



T < T E K (t) + Ce < T e (^)(0) + Ce 



< 



< 



and an analogous inequality holds for v~. Therefore vt(t), v^(t) G Here 
we denoted with T £ the sum 1/ + CoQ defined on e-approximate wave front 
tracking solutions (see [H formulae (7.53), (7.54)]). We may thus apply 
Proposition 13. 12l Proposition 13.141 and the main result in [U Chapter 8], 
that is [H Theorem 8.2], to obtain 

S(v £ v (t),v £ v (t)) = ^(v £ v (t),v £ (t) 



v £ , vf^j (0) + Cet + Ce v e ri {t) - v £ v (t) 



< §(v v , v n ) + Cet + Ce v^t) - v £ (t) 



L 1 



+ Ce\\v v - v v \\ L1 . 
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Recall that as e — > by [4] Theorem 8.1] v^{t) — > Stv v and v^{t) — > 5(W^. 
Hence, Proposition 13.131 and (|3.8p ensure that 

"3(S t v v ,S t v n ) < liminfS ^(i),^(t)J < * («,,,«,,) < S(u,tt) + 77. 

By the choice of and we have that v v — > u and v v —> u in L 1 as 
7/ — > 0+. Therefore, using the continuity of the SRS in L 1 and applying 
again Proposition 13.131 we may conclude that 

S(S t u,S t u) <limmiS(Stv v ,S t v v ) < H(u, u) , 

T)— >0+ 

proving f££j. The latter item fiii,) follows from Proposition 13.131 □ 



4 Wave Measures Formulation 

Let / satisfy (F) and u € T>$ as defined in (|2.6|) . Since TV(u) is bounded, by 
possibly changing the values of u at countably many points, we can assume 
that u is right continuous. Its distributional derivative /i is then a vector 
measure that can be decomposed into a continuous part /j c and an atomic 
one fi a . For i = 1, . . . , n, consider now the wave measure 

Vi{B)= I k{u)diJL c + ^E i [u(x-),u{x+)) (4.1) 

where li(u) is the left i-th eigenvector of Df(u), Ei is the i-th component of 
the map E defined at (|2.3|) and B C R is any Borel set. Here and in what 
follows, we assume that l\, . . . , l n are normalized so that 



k(u) ■ Tj(u) 



1 if f = j 
ifi ± j 



with rj as in (|2.ip . Let respectively /j~ be the positive, respectively 
negative, part of /ij and = fif + fi~ be the total variation of /ij. Aiming 
at a definition of the interaction potential by means of the wave measures, 
introduce the measure 



p= Yl Ia**i ® |a*j| + E(^< +pi®Pi ®^ + ) ( 4 - 2 ) 

l<j'<i<n i=l 

and, as in (3J HI [5], set 

Q(«) = p({(z,y) £l 2 :x<i,}J (4.3) 

n 

T(u) = 53|/ii|(R) + Co.Q(«). (4.4) 



i=l 
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For u, u in T>s, we now define the functional 

n ^+00 



" p+00 

S(u,u) = ^2 \qi(x)\Wi(x)dx 
i=i J -°° 



(4.5) 



where the weights Wj are defined by 

Wi(x) = 1 + KiAj(x) + Ki k 2 (Q(u) + Q(u)) . (4.6) 

Here, n\ and K2 are as in [H Chapter 8], see also (|3.2|) . By means of the 
wave measures \ii and /ij of u and u, if the z-th field is linearly degenerate, 
define the weights Aj by 

A i0) = M O -00 '^) + (]-«)>«]) 

l<i<j<n 

+ X] W (] x ' +oo D + 1 A? I (R+°°Q 

l<j'<i<n 

whereas in the genuinely nonlinear case we let 
A* 0*0 = M _oo ' x ]) + (]-°°^]) 

l<i<j<n 

+ W a; '+ oo D + (] x ' +oo D 

l<j<i<n 
( n 

(N {}-oo,x}) + (]x,+oo[)\ ifqi(x) < 



+ < 



i=l 
n 



(]x,+oo[) + \fj,i\ (]-oo,x])j if qi(x) > 



i=i 

On Dg, Q and T obviously coincide respectively with Q and T, therefore 
also H, <1? and S all coincide on T>%. Below, we prove that T = T, Q = Q 
and S = S on all T>s 

The following result is a strengthened version of [H Lemma 10.1]. 

Lemma 4.1 There exists a positive C such that for all u G "Dg, i = 1, • ■ ■ ,n 
and a, b £ M with a < b 



Ei (u(a+),u(b-)) - m (]a, b[ 



< C ■ diam 



u[]a,b[)) ■ H (KM) 



Proof. We use below the following estimate, see |H p. 201], valid for all 
states u, u, 

\Ei(u, u) — k(u) ■ (u — u)\ < C \u — -£t 1 2 . (4-7) 
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By the triangle inequality, 



Ei (u(a+),u(b-)) (]a,b[) 



< 



Ei («(a+),u(6-)) - k{a+) ■ (u{b-) - u(a+)) 



+ 



Hi (]a, b[) - k(a+) ■ (u(b-) - u(a+)) 



The first term in the right hand side is bounded by (|4.7p . By (|4.ip . for 
I = ]a,b[. 

(M(I) = J (li (u(0) - k («(<*+))) • dfi c (0 

+ £ (Ei («(£-), «(£+)) - «i («(£-)) ' («(£+) - 

+ £ (*« ( n ^-)) - • («(£+) - u (t-)) 

+k (u(a+)) • ^£ («(£+) " «(£-)) + / ^c(0 j • 

We now estimates the different summands above separately. The Lips- 
chitzeanity of k ensures that the first summand above is bounded by 



/ (k («(£)) - k («(«+))) • rfMc(0 < C • diam («(/)) • . 
Passing to the second summand, using (|4,7|) 

£ (Ei (u(£-),u(£+)) ~ U («(£-)) ' («(£+) " "(£-) 

< c£|u(£+)-u(£-)| 2 

< C • diam («(!)) • |//| (/) . 

Using again the Lipschitzeanity of Zj, the third summand is estimated as 



£ (ii Me-)) - ^ («(«+))) • («(£+) - «(*-)) 

while the last one can be rewritten as 



< C ■ diam («(/)) • 



(«(a+)) • ( £ («(£+) - «(£-)) + ^ dMc(0 
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k (u(a+)) ■ [n a (I) + n c {I)) 

k («(«+)) • n(I) 

U (u(a+)) ■ (u(ft-) - u(a+)) 



completing the proof. 



□ 



Lemma 4.2 Let u £ T>$ with wave measure fjL\, . . . Then, there exists 
a sequence v u € T>* & with wave measures /i", . . . , n v n such that 



lim \\v u — u\\ L c 



0. 



lim \\v u — u\\ T i=0 



lim ^ (J) = fJ-f(I) for any interval I C 



lim Q(v u ) = Q(u) . 

v— >+oo 

Moreover, an explicit definition of such a sequence is (|4.10p . 
Proof. For z/ 6 N \ {0} choose a, b € R with 6 — a > 1 and so that 



(4.8) 
(4.9) 



I A* | (]-oo,a]) + |/x| (|>,+oo[) < - 

r+oo 

\\u(x)\\dx+ / ||tx(x)|| dx < — 



Choose a finite sequence of real numbers x\, . . . , xn so that a = x\ < x<i < 
■ ■ ■ < xn-i < xn = b and (]x a -i, x a {) < / b _ a \ v ■ Introduce the points 
Vo = xi- 1, y a = (x a + x a+ i)/2 for a = 1, . . . , N — 1 and vn = xn + 1- Let 

v v = £ («(*«-) X [ya _ li2a[ + «(*«+) X K)fe[ ) • (4-10) 

Due to the above definitions, the L 1 and L°° convergence v u — > u is 
immediate (observe that both v v and u are right continuous). 

We first consider the intervals / = ]— oo, x] for x 6 R or I = R. Let /x^ 
be the wave measure corresponding to v u . For notational simplicity, below 
we set xo = — oo, x^ + \ = +oo, u{xq+) = and u{xn+\— ) = 0. Let a be 
such that x G [x^, (a = N for / = R). For any i = 1, 



,n, 



l/^ICO - | Mil CO = ^ (tt(x a -),u(x a +)) 

+ ^2 \Ei (u(x a -i+),u(x Q -)) 



(4.11) 



a=l 



+ \Ei (u(Xa+),u(Xa+l-)) 
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Q = l 



Observe that 



Ei (u(x a -),u(x a +)) = ^ ({x a }) for all a = l,...,N 

Ei (u(x a -),u(x a +)) = \fii\({x a }) for all a = l,...,N 

Ei (u(Xa + ), u(Xa+l-)) < C \fl\ (}Xa, Xa+l[) < C/v 

\Hi\(]Xa,x]) < C\fJ,\ (\Xa,X & +l[) < C/v 

so that, by Lemma 14.11 



ft 



C 
v 



|(J) - \ih\{T) = \Ei(u(x a -i+),u(x a -)) - \fJ>i\ (]x a -i,X a [) j + 

a=l ^ ' 

< y](\Ei(u(xg-i+),u(x a -)) - m (]x a -i,x a [) J H 

a=l ^ J v 



< ^ \Ei (u(x a -i+),u(x a -)) - m (]x a -i,x a [) 



a=l 



c 
+ — 

v 



< C ■ y~]diam [u [\x a -i,x a 

a=l 

51 1 c 

< c-y^-lnl (]x a -i,x a [) h — 



Xct—l 



,x a [) + 



C 



< C.(l + H(R)).i 
and limsup uuf (J) < \Hi\{I)- 

v— >+oo 

Passing to the other inequality, introduce the functions 

Wi(x) = Hi (] — °°> ^]) an d tf^(x) = /ij' (]— oo, x]) . 
and repeat the same computations used in (|4.11|) - (|4.12|) to obtain 



(4.12) 



\x) - Wi (x)\ = \hUi) - IH(I)\ < C ■ (1 + H W) 



1 



showing that — > Wi uniformly on R. By the lower semicontinuity of the 
total variation 

=TV (wi,]-oo,x\) <liminfTV «,]-oo,x]) = lim inf I /tf I (/) 
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showing that Lu^ (/) — > |/Uj|(7) as v — > +00. This convergence, together 
with the uniform convergence above, imply that 

lim U> + (I)+rt-{I)) = 

v— >+oo V / 

lim = /x+(I)- M r(J) 

which together imply (|4.8p for / = ]— 00, x] or / = 1BL Let now u(x) and v v {x) 
be the right continuous representative of respectively u{— x) and v u (— x) with 
corresponding wave measures /2j and ft". The previous computation show 
that (]— 00, x]) — > /i^ 1 (]— 00, x]) for all real x. But (]— 00, x]) = 

IH' ([ — X )+°°D an d fit (]— 00, x]) = /U^ ([— x, +00 [) . Therefore, (14, 8p holds 
also for the intervals / = [x, +00 [ and therefore for all real intervals ICR. 

Passing to (14. 9|) we observe that it is enough to show the convergence 
of every single term in the sum (|4.2p which defines the measure p. Since 
the computations for these terms are identical, we show the convergence of 
only one, say (g> p"'~ ■ Fix e > and choose a finite set of real numbers 
xi, . . . , xpf e so that — 00 = Xq < x\ < ■ • • < x £ Ni: < x £ Ns+l = +00 and 

I /la I ^q.!,^ [j < e. To simplify the notations, we define 

K = I (x, y) G R 2 : x < y} , r 1 ' = <g> pf and r = ® ftr. 
Now, write -ftT as the union of a finite family of disjoint sets as: 

Ra,f3 = [ x a' x a+l [ X 

x a, /3=o ' a=0 

and compute 

N e 

\t»{K)-t(K)\ < Yl \r v {Ra,p)-r{R^)\ 

a, = 
a<f3 

+ ^ {r u (R a ,a DK)+ T(R a , a n K)) . 

a=Q 

We have also the estimate 

T V (R a>a KK) < T»(R a , a \{(x £ al x%)}) = T»{R a , a )-T»({(x%,x e a )}) 

r(R a , a nK) < T(R a , a \{(x%,x%)}^ = T(i? Qia )-r({(x^,x^)}) 

The limit (|4.8p implies that the product measure converges on rectangles: 

lim T u (R a p) = t (R a ,p) 



h fii x /3+l 



K = ( (J i? Ql/3 U |J (i?a,a n 
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Therefore 



\imsup\r^K)-r(K)\ < 2 £ ( r{R a , a ) - r ({(x £ a , x E a )})) 



Q = 



a=0 



2E K Kl ^(Ka.<+iD 



+^ + (]<><+lQ ^ ({<}) 

< C^eH ([a£,a£+i[) 

< Ce|//|(R). 

By the arbitrariness of e, lim t v (K) = t(K), completing the proof. □ 

v— >+oo 

The following proof, simpler than that in [4, Theorem 10.1], is based on 
the piecewise constant approximations introduced in Lemma l4.2i 



Proposition 4.3 The Junctionals Q and T defined in (|4.3p and (|4.4p are 

Zower semicontinuous with respect to the L 1 norm. 

Proof. We prove the lower semicontinuity only of T since the other one is 
easier. Take u v , u G T>s such that u v — > u in L 1 . By Lemma 14.21 there 
exists a sequence of piecewise constant functions v u € V* & such that 

IK _ ^IIl 1 - ~ ' ^i v v) — "^(^) + ^ • 

If we define Z = liminf^ +00 T(u^), then Z < lim inf j,_> +00 "Y(u v ). Moreover, 
by possibly passing to a subsequence, we suppose that Z = lim^+oo "T(v u ) 
and that v v {x) — ► u(x) for every x € D, where M \ D has zero Lebesgue 
measure. Fix now e > and, using again Lemma 14.21 choose a piecewise 
constant function v £ € D| which approximates u: 



N 



^ ( \ 

E\ u(x% — )y, r + w(a£, + ) Y r r) 

a=l x ' 



r(u) <T(v £ ) + e = T(v £ )+e. 
By Remark 12. 1| for every a = 1, . . . , N e , we can choose points 

Ce]ti.4[nfl and d%+ e]a£,i£[nz> 
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such that the function 

a=l x ' 

satisfies T(v £ ) < T (v £ ) + e. Define now 

Since — > n fe, ), we can apply again Remark 12,11 to obtain that 
for v sufficiently large one has T (v £ ) < T (v £il/ ) + e. 

But v £ ^ is obtained by removing an ordered sequence of values attained 
by v u , therefore we can apply Proposition 13.21 to get T (v £tU ) < T (v u ) and 
therefore we have the following chain of inequalities: 

f (u) < T (v E ) + e < T (v £ ) + 2e<T (v Ejl/ ) +3e<T (y v ) + 3e . 

Taking the limit as u — > +00 one obtains T (u) < l+3e and the arbitrariness 
of e > implies T (u) < I < lim inf 2,^+00 Y(u u ). □ 



Corollary 4.4 T/ie junctionals Q and T defined by (I4.3P and (|4.4p coincide 
on all T>s with Q and T . 

Proof. It is a straightforward consequence of the lower semicontinuity of 
both Q and T (Proposition I4.3P and Q and T (Proposition 13. 1 j) . 

Indeed, consider only T and T. They obviously coincide on Fix now 
u € T>s- By the definition (|3.ip of T, there exists a sequence U;, of functions 
in £>| converging to u in L 1 and such that T(u^) — > T(u) as ^ — > +00. By 
the lower semicontinuity of T (Proposition I4.3P , we obtain 

T(n) < lim inf Y(v u ) = lim inf Y(v u ) = lim T(v u ) = T(u) . 

v— H-00 >+oo z/^+oo 

Analogously, by Lemma 14.21 we can take a sequence v u of functions in T>% 
such that v v — > n in L 1 and T(iv) — > T(n) as f — > +00. Therefore, along 
this particular sequence, we may repeat the estimates as above applying the 
lower semicontinuity of T (Proposition 13. IP : 

Y(u) < liminf Y(v v ) = liminf T(v u ) = lim T(v u ) = Y(u). 

V— > + 0O i^^ + OO V— > + 00 

□ 
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By Corollary 14.41 m the following we write Q and T for Q and T. 

The following proposition shows the lower semicontinuity of S along 
piecewise constant converging sequences. 

Lemma 4.5 Let u, u G T>$. Then, the approximating sequences v u , v u E 
V* & defined in Lemma \JJ^ satisfy also lim^ +00 S(v 

Proof. Define q(:c) and cf{x) so that u(x) = S (q(x)) (u(x)) and v u (x) = 
S (q i/ (x)) (v v (x)). Then, q 1 ' — > q uniformly and in L . Let Wj(x), re- 
spectively W^(x), be the weights defined in (|4.6I) with reference to u, u, 
respectively^,^. Compute 



\v v ,v v ) 



[u,u, 



< 



E 

8=1 
n 

+E 

i=l 



\qi(x) — qi(x)\ -W^(x)dx 

oo 

\ qi (x)\ ■ \ W^(x)-W i {x)\dx 



The first integral converges obviously to zero. Concerning the second one, 
where qi{x) = the integrand vanishes. Otherwise, if qi{x) ^ 0, then for 
v sufficiently large qi(x) ■ q%(x) > and hence the weights depend continu- 
ously only on the wave measures and on the interaction potentials which all 
converge by Lemma 14.21 Therefore, for all x € M the integrand satisfies 



lim \q i (x)\-\W^(x)-W i (x)\=0. 



The Dominated Convergence Theorem concludes the proof. 



□ 



Proposition 4.6 For any u, u £ T>$ and v u , v u € 2?| such that v v — > u 
and v v — > u in L 1 as v — > +oo, we have S(u,u) < lim inf ,,^ +00 &(v u , v v ). 

Proof. Let I = liminf^ +00 S(v u , v u ). Passing to subsequences, we assume 
that I = lim^+oo H (v u , v u ) and that v u , respectively v u , converges pointwise 
to u, respectively u, on a set DC1 with M\D having zero Lebesgue measure. 
Fix e > and apply lemmas 14.21 and 14.51 to find two functions 

Ne ( \ 
v e = J2 {<<-)x [yi _ iiX%[ + u(xl+) X j , 

a=l x ' 
a=l x ' 
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such that H(u, u) — S(v e , v £ ) + \\u — w e || L i + \\u — We|| L i < e. On D|, S 
and $ coincide, hence Remark 13.51 applies and there are points 

x £ a &]y £ a -i,x £ a [nD and x e + e]x £ a ,y £ a [nD 
Cejiia-D^tnf and ^ + 6]i £ a ,!/a[nD 
such that the two functions 



a=l 

£ 

a=l 



satisfy S(v e , € e ) - S(# £)i „ « £jW ) + ||u e ^ - u e || Ll + 
sufficiently large. Lemma 13.91 thus implies 



< e for v 



K^e,v,v e , v ) < a(v u ,v u ) + C \\v £ ,u ~ v v L, + 



Hence 



E(u,u) < E(t> e ,u £ )+£ 

< S(-y e>I/ , v e)U ) + 2e 

< 3{v u ,v u ) +2e + C (\\v £ ,u - VuW^i + 

< £(«„,«„) + 2e ■ (C + 1) + C (\\v e - v„|| L i + ||u 6 - u„|| Ll ) 
so that, for v — ► +oo, 

S(u, u) < Z + 2e(C + 1) + C (\\v £ - it|| L1 + ||u e - u|| L i) < / + 2e(2C + 1) . 
The arbitrariness of e concludes the proof. □ 



Theorem 4.7 Let f satisfy (F). The functional a defined in (|4.5p coincides 
on all T>s with S as defined in (|3.4|) . In particular S is lower semicontinuous. 

Proof. Both functionals can be approximated through their evaluation on 
piecewise constant functions (see (I3.4D and Lemma 14.51) . Both function- 
als coincide on piecewise constant functions and are lower semicontinuous 
along sequences of piecewise constant functions (see Theorem l3.6l and Propo- 
sition 14. 6|) . A procedure identical to that of Corollary 14.41 completes the 
proof. □ 
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